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SOME ENERGY CONSIDERATIONS IN GAMMA RAY BURST LOCATION 
DETERMINATIONS BY AN ANISOTROPIC ARRAY OF DETECTORS 


BY 


John H. Young 

Associate Professor of Physics 
University of Alabama at Birmingham 
Birmingham, Alabama 


ABSTRACT 

The anisotropic array of detectors to be used in the Rurst 
and Transient Source Experiment (BATSE) for locating gamma ray 
burst sources is examined with respect to its ability to locate 
those sources by means of the relative responses of. its eight 
detectors. It has been shown that the energy-dependent 
attenuation effects of the aluminum window covering each detector 
has a significant effect on source location determinations. 
Location formulas have been derived as a function of detector 
counts and gamma ray energies in the range 50 - 150 keV. 

Deviation formulas have also been derived and serve to indicate 
the location error that would be caused by ignoring the influence 
of the passive absorber. 
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INTRODUCTION 


One of the four experiment packages aboard the Gamma Ray 
Observatory (GRO) , scheduled for release into orbit from a 
shuttle flight in May 1988, is the Burst and Transient Source 
Experiment (BATSE). Utilizing an array of eight Nal(Tl) 
detectors BATSE will, as one of its major objectives, attempt to 
determine the angular locations of gamma ray burst sources by 
means of analyzing the relative count responses of the detector 
array. Figure 1 depicts the GRO configuration^. If the planes 
of the detector faces, each having an area of about 2000 cm , 
were extended so as to form a closed surface, that surface would 
be the regular octohedron also depicted in Figure 1. 

The 1.27 cm thick detector crystals will be overlaid by 2mm 
of aluminum. The gamma ray flux incident on the aluminum windows 
will be attenuated by that material and a diminished flux will 
consequently be incident on the detectors. The relative 
orientation of the detectors, dictated by the octohedral 
geometry, necessitates that gamma rays incident on one detector 
will have transversed a different thickness of aluminum than 
those incident on that detector's neighbors. Thus, in addition 
to being energy dependent due to the material's attenuation 
coefficient, the attenuation will be orientation dependent. 

The work described here is the result of an effort to take 
the energy dependence of the relative detector responses to a 
burst into account in the angular location determination of the 
source. The energy-dependent attenuation effects are discussed 
and location formulas containing these effects are presented. 
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THE BATSE GEOMETRY 


The BATSE detector array conforms to the qeometry of a 
regular octohedron shown in Figure 2. The planes, corresponding 
to the eight detector faces, are identified by the numbering 
scheme indicated. The outward normals from the eight planes are 

A A A 

given by the following, in which i, j, and k are unit vectors 
along the positive x, y, and z axes, respectively: 


A 1 A A A 

U, = ±(i + j + k) 

-/3 


A 1 A A A 

u 2 = i (i + j " k > 

J3 


A A 


= -(-i + j + k) 

v3 


A A 


u 4 ’ jt' 1 + J " k) 


Ur = p( -i - j + k) 
D v3 


Ur = 




- j - k) 


U 7 = 


1 A A A 

“(i - j + k) 

n 


A I A A A 

u R = i - j - k) 

8 s 


in 


The following relationships between these normals are worthy of 
note and will be used later: 


A AA A A A A A 



A A A A 

u 4 = u 3 + u 2 “ U 1 * 


If the octohedron is imagined to be centered at the origin 
of a large sphere of radius R (Figure 3) from whose surface a 
burst originates with angular coordinates 9 and <j>, the detected 
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portion of the burst flux will be that which propagates along the 
direction -$/R. Let t represent the flux density (photons per 
unit area) in this direction, with magnitude I Q ; it can be simply 
represented as 

. ^ A A A 

I = I (- -=r) = - I (sin8 cosd>i+sin0 sin<b j + cos0k) . (3) 

OR o 

BATSE data of interest here are derived from the incidence 
of f onto any combination of the eight planes of the 
octohedron. The number of photons passing through the i-th plane 
whose area is A will be 

N. = I •(- u.A) = I A cosa. , (4) 

1 1 O 1 

A 

where a. is the angle between the burst direction and - u. . 

l a l 

DETECTOR RESPONSE 

Only under ideal, and thus rarely attainable, conditions 

A 

would a detector whose outward normal is u^ actually record the 
number of photons given by Eqn. (4). That expression must be 
multiplied by the detector efficiency, resulting in a decrease 
from Eqn. (4). The efficiency is an energy-dependent factor 
which accounts for attenuation losses in the passive absorber 
covering the detector crystal as well as the probability that a 
photon will interact with the crystal upon entering it. If t and 
D represent, respectively, the window and crystal thicknesses, 
the efficiency for the case of normal incidence is given by 
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where y A (E) and y^(E) are the energy-dependent attenuation 
coefficients of the aluminum and sodium iodide. If, however, a 
photon were to enter the system at some anqle a relative to the 
system axis, it would transverse a thickness t/cosa of absorber 
before getting into the crystal. Further, propagation at that 
angle gives the photon an increased effective thickness of 
crystal with which to interact and thus be recorded. The latter 
effect amounts to replacing D by D/cosa in Eqn. (5). The 
detector efficiency will be used here in the form 

-y. ( E) t/cosa -y LE ) D/cosa 

n(E) = e A (1 - e ) . (6) 


The number of photons, , recorded by the i-th detector 
will then be given by the product of Egns. (4) and (6); that is 

-y (E)t/cosa. -y (E)D/cosa. 

NN (E , 0 , <f>) = I o Acosa.e (l - e ) . (7) 

The energy dependence of the considered here will be the 
result of the attenuation coefficients' dependence on photon 
energy. The relationship of to the angular location of the 
burst is dictated by the functional relationship between a^ and 
( 6 , <f> ) ; these relationships are easily determined from Eqns. (1), 
(3), and (4). The objective of the present work is to develop 
the means of determining 0 and <j> in terms of the recorded 
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It is apparent from Eqn. (7) that the 


detector counts, N^. 
computation of the angular coordinates of the burst will give 
results dependent on photon energy. The following discussion 
will, in particular, be focussed on the sort of corrections 
required in location determinations mandated by the inclusion of 
energy considerations above 50 keV. 

“U t / COSQl 

The exponential function e has maximum value 

(0<a<90°)e ^ , corresponding to normal incidence. Table 1 lists 
the attenuation coefficient values for aluminum over the energy 
range of interest (50 - 200 keV) and compares the value of the 
exponential e for t = 0.2 cm with its linear approximation, 

1 - pt. The largest discrepancy between the exact and 
approximate values is seen to occur at 50 keV, at which the 
difference is about 2%. At 200 keV, on the other hand, the 
linear approximation of the exponential function is seen to be 
accurate to within 0.2%. It will be shown below that the most 
probable energy of a photon within a typical burst is about 100 
keV. The linear approximation agrees with the exponential at 
that energy to with 0.5%. It would appear reasonable then to 
approximate the first exponential appearing in Eqn. (7) by its 
linear representation, and thereby obtain angular locations which 
will be "first-order," a refinement over the "zero-th order" 
attainment based on totally ignoring attenuation effects. 

Table II displays values of attenuation coefficients of 
sodium iodide^ over the energy range of interest, along with 
values of the second exponential factor appearing in Eqn. (7). 

It is seen that for a 1.27 cm-thick crystal, this factor is 
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essentially unity up to about 125 keV. It is important to 

recognize, however, that the values displayed in the fourth 

column correspond to normal incidence and thus correspond to 

minimum values of (l - e ' jD / cosot ^ > Alternatively, that column 

can be looked upon as exhibiting the proximity of this factor to 

unity. More typically an angle of incidence a ** 45° might be 

encountered, thus requiring photons to traverse a greater 

thickness of the crystal, while increasing their probability of 

detection. The latter values are tabulated in the fifth 

column. It should be noted here that if these same 

considerations regarding the angle of incidence were given to 

Table I that the exponential factor and its linear approximation 

* 

would be in agreement to within 0.8%. The burst location 
discussion presented here will assume validity of the linear 
approximation of attenuation in the aluminum window and that the 
crystal will respond to any photon which enters it. 

Gamma rays within the 50 - 150 keV then can be expected to 
give rise to the following approximate forms of Egn. (7): 

N.(E),e,*) * I q A( cosa. - u A <E)t) (8) 

The fact that the attenuation coefficient of sodium iodide no 
longer appears in this approximation allows the subscript to be 
dropped from hereafter p(E) will refer to aluminum. 
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BURST LOCATION! DETERMINATION! 


Use of Eqn. (8) in determining the angular coordinates of a 
burst source will now be described. The energy dependence of the 
results will be emphasized and location formulas including the 
energy dependent term in Eqn. (8) will be derived. 

Consider a burst activiting detectors 1 through 4. (This 
means that 55° < 9 < 125° and 45° < 4> < 135°. Generalization of 
the results to other detector combinations will be discussed 
later.) Using the definition of cosa,. = ( ?/I q ) - ( — u ^ ) along with 
Eqns. (1) and (3) gives the following expressions for 
the cosa^ in terms of 9 and <J> : 

cosa. = ^(sin9 cos(i + sin0 sin$ + cos0) 

V3 

cosa 2 = yj(-sin0 cos$ + sin9 sin<|> - cos0) , 

1 (9) 

cosa- = — (-sin0 cos$ + sin0 sin<j> + cos0) , 

V3 

cosa. = — (-sin0 cos<f> + sin0 sin<J> - cos0) . 

V? 

At this point two important identities between the cosa. 
should be noted. It can be immediately seen from inspection of 
the statements given in Eqn. (9) that 

cosa^ = cosa^ + cosa 2 ~ cosa^ . (10) 
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This relationship simply restates the last connection between the 

a 

ik noted in Eqn. ( 2 ) and indicates the sufficiency for the use of 
three detectors only. A less obvious but, as will be seen below, 
very important additional identify which can be easily verified 
is the following: 

2 4 2 2 

(cosa^ - cosctjj = -k - [(cosc^ + cosa^; + (cosa^ - cosa^) J .(11) 

It is reasonable to expect that burst location determinations 
could be made with increased confidence for the case of the more 
highly activated detectors. Eqn. ( 10 ) will allow those 
determinations to be made in terms of the three highest detector 
responses, assumed here to be N-^, N2, and N3. The humber of 
photons incident from a given burst, IqA, appears in each 
detector statement of the form given in Eqn. (8). This cuantity 
is unknown, however, and Eqn. ( 11 ) serves the purpose of allowing 
it to be determined. 

Combining Eqns. ( 9 ) and (8) gives the following system: 

= il o A[sin8 cost + sin 9 sin<)> + cos 9 - y(E)tl 

= il A(sin 9 cost + sin8 sin<}> - cos 9 - y(E)t] ( 12 ) 

V3 0 

N_ = —I A[-sin8 cos<t + sin8 sin<t + cos 9 - y(E)t) 

J V3 0 

The explicit statement for will be suppressed since Eqn. ( 10 ) 
implies that N 4 = N3 + N2 - N^. It is more convenient to express 
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the angular functions of Eqn. (12) in terms of a "corrected 
count" expressed as 


Then 


N. ’ = N . + — I Au(E)t ; 

1 1 3 o 


N, ' = — I A(sin0 cosiji + sin0 sin<{> + cos0) , 

1 75 ° 


( 13 ) 


N ’ = — I A(sin0 cos<j> + sin0 sin<|> - cos0) , 

* 75 ° 


(14) 


N 3 1 = — I Q A(-sin0 cos4> + sin0 sin<{> + cos0) 


Eqns. (14) are easily combined to yield expressions which 
can be solved for 0 and <t> . It is apparent that 


V + V 


— I A sin© sin<(> , 

V5 0 


and 


(15) 


V ~ V 


— I A sin0 cos 6 . 

/5 0 


(16) 


Division of these and use of Eqn. (13) gives 


tan<j> 


w 2 * N 3 
N ‘l * N 3 


N 2 + N 3 ,VJ r o Au(E)t 
N 1 - N 3 N 1 ' N 3 


(17) 


The second term on the right of Eqn. (17) will be viewed as a 
connection to the azimuthal location of the burst, and is due to 
the finite thickness of the detectors' aluminum window. Its 
significance will later be discussed in detail. 
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The expressions for determining 0 and <{> in terms of N 2 , 

and N 3 are seen to be energy dependent and contain the unknown 
flux density I Q . The latter quantity can be eliminated from 
Eqns. (17) and (21) by expressing it in terms of the detector 
counts. This can be achieved by use of the cosa. identity given 
in Eqn. (11). In terms of the N^' f 


N ' N ' 2 . N ' N ' 2 N ' N ' 2 

~ Y~K ) 3 ^T"A + I~A^ + ^T~A “ I ^ 

00 00 00 


( 22 ) 


If the N^' are expressed according to Eqn. (13) a quadratic in 
I q A results which has solution 
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C - II 



2u(E)t(N 2 + N 3 ) 


I A = 


[4 - 4(p(E)t) 2 ] 


{[2y(E)t(N 2 +N 3 )] 2 + [4 " 4 ( y ( E) t > 2 ] [(t^- + N 3 ) 2 + (N 2 +N 3 ) 2 ] } 


[4 - 4(p(E)tr] 


(23) 


The negative root has been discarded by the requirement I q A > 0. 

Equations (17) and (21), along with (23), are sufficient to 
determine the angular location of the burst source for photons of 
energy between 50 keV and, perhaps, 150 keV. The expressions to 
which these quantities reduce in the limit t -*■ 0 are those which 
can easily be derived by ignoring attenuation (and .thus energy) 
effects from the outset. The latter expressions, as previously 
mentioned, can be thought of as "zero-th order expressions" and 
will hereafter be superscripted "(o)". Explicitly, 


tnnfi (o) - N 2 + N 3 
tan<}> - _ N t 

1 3 


(24) 


tan0 


o) ^ N l 2 + N 2 2 + 2N 3 2 + 2N 3 (N 2 “ N l^ 1/2 


N 1 - N 2 


(25) 


and 


‘V 0 


(o) _ 


A t (N i 


N 2 ) 2 + (N x - N 3 ) 2 + (N 2 + N 3 ) 2 ] 


^2 


( 26 ) 


That neither Eqn. (24) nor (25) contains the flux density 
explicitly means that Eqn. (26) need not be evaluated for 
locations in this order. 
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The effect of the aluminum window on burst location 
determination is seen, by coraparinq Eqns. (17) with (24) and (21) 
with (25), to give rise to an abberation. The role of the energy 
dependence in this abberation will now be examined. Computations 
of the source location and flux density must be made in terms of 
detector counts. The energy influence on these quantities can 
perhaps best be illustrated by choosing reasonable values of what 
might typify data registered by the three detectors. For this 
purpose let = 3, Nj = 2, and Nj = 1. (These numbers have, of 

3 

course, been scaled down by '10 , which is permissible since only 
the relative counts really matter.) These values give, with t = 
0.2 cm and Eqns. (23), (17), and (21) the results displayed in 
Table III. The deviations, 6 4> = <J> - and 50 = 0 - 0^°^, are 

seen to diminish with increasing energy. This result is what 
would be expected, as the attenuation would be reduced for the 
more energetic photons. It would appear, however, that in no 
case displayed there is the deviation negligably small. 

Although not small with respect to burst location 
determinations 6$ and 60 are small in the sense that their 
tangents can be set equal to these angles (in radians) 
themselves; e.g., tanS$ = 5 <J> . These approximations are known to 
hold for angles considerably larger than those displayed in Table 
III and allow general deviation formulas to be derived. 

According to Eqn. (17) 


tan<)> - tan<|> 


(o) V? r o A * 1E)t 
N 1 - N 3 
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Setting <J> = + 6<J>, using the addition formula to 

express tan(<J>^°^ + 6 4> ) , with tan64 » 6 4> gives 


tanift^ 0 ^ + 64> 
I- 64 . tan<(> (o) 


tan<|> 


(o) A 1 o Aw(E)t 
N 1 - N 3 


Solving for 64 and expressing tan<|>^ 0 ^ according to Egn. (24) 
gives the result 


6 <f> = 


-7T 1 o 


I_Ap(E)t(N 1 -N 2 ) 


(N 1 -N 3 ) 2 + (N 2 +N 3 ) 2 +j= I Q Au(E)t(N 2 +N 3 ) 


(27) 


As might be anticipated from Eqn. (21) the corresponding 
calculation of 66 is more involved. The fact that tan0 is 
expressed in terms of a square root necessitates solving a 
quadratic equation in 60. The somewhat laborous details of the 
calculation will be omitted here in favor of getting on with the 
result : 


60 = ~ tan9 (o) [l+C+tan 2 0 (o) ] 

1 [1 - Ctan 2 0 (o) - tan 4 0 (o) ] 


{tan 2 0 (o) [l+c+tan 2 0 (o) 1 + C [ 1-Ctan 2 0 (o) - tan 4 0 (o) ]} /2 

[1 - Ctan 2 0 (o) - tan 4 0 (o) ] 


where tan0^°^ is given by Eqn. (25) and 


(28) 


C = 


A I Q Au(E)t(N 2 +N 3 I q Au (E)t) 

(Nj-Nj ) 2 


( 29) 
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The choice of sign is made according to the requirement that 
60 ± 0 as t 0 (i.e., as C + 0). This means that if 

tan0^°^ > 0, choose 50 + ; if tan0^°^ < 0, choose 60_. It might 
appear that 60 would become infinite for that value of 0^°^ which 
causes the denominator to vanish. It must be kept in mind that 
the above was derived under the premise of the activation of 
detectors 1 through 4. This means that 55° < 0 < 125°; it is 
easily shown that the value of 0^°^ which would cause the 
denominator of Eqn. (28) to vanish is not within this range. 

The fact that the deviation formulas are energy dependent 
makes their expressions, Eqns. (27) and (28), somewhat difficult 
to interpret. (It should be kept in mind that the appearinq 
in those expressions are themselves effected by gamma ray energy: 
the fact that p decreases with increasing energy means, according 
to the relations given in Eqn. (12), that the will correspond- 

ingly increase and thereby additionally influence 6 <|> and 5 0.) 

The significance of these formulas might best be appreciated by 
expressing the deviations in terms of the most probable energy of 
a photon in a burst over the 50 - 150 keV range. 

Burst spectra so far obtained appear to obey the following 
empirical relation^: 


dn 1 
dE “ e 1.4 * 


The most probable energy, <E>, is then calculated as 
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<E> 


J n ( E ) EdE 
j n ( E ) dE 


t 


(30) 


integrated over the afore-mentioned limits. The result for the 
above power law gives <E> = 96.4 keV. 

The deviations given by Eqns. (27) and (28) can now be 
assessed at this value of energy for various detector counts. 
Table IV displays a variety of detector response examples along 
with the corresponding deviations at 96.4 keV. (The value 
of u used here was 0.513 cm - ^.) The column headed n (66 ± )" has 
been included to indicate the choice of sign used in Eqn. (28). 
The zero values in and 66 result from the equality of Ng and 
Ng, and Ng and N 2 , respectively. 


DISCUSSION AND CONCLUSIONS 

Burst location determinations based on detector responses 
have been discussed in a way that includes gamma ray energies. 
Equations (17) and (21), along with (23) can be used for energies 
in the approximate range 50 - 150 keV. Although those 
expressions were developed under the assumption of a burst 
activating detectors 1 through 4, so that at least two of these 
have N g > 0, it is a straight-forward matter to qenerate the 
appropriate expressions in the event of activation of some other 
detector combination. The unit vector relationships given in 
Eqn. (2) indicate, e.g., that N 7 = - N^ and Ng = - Ng. Thus, if 
a burst activated detectors 1, 2, 7, and 8 it would only be 
necessary to replace Ng in the expressions given by -Ng. N 2 need 
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not appear in the final expression, as was the case with in 
the situation discussed, since N 2 = N 1 + N 8 “ N 7* The deviation 
formulas can be similarly altered. 

The examples of deviations associated with the detector 
responses considered in Table IV assumed a particular power-law 
energy spectrum of the burst in order to compute the most 
probable energy, <E>. This quantity can easily be calculated 
from Eqn. (30) for any burst spectrum. If, for example, the 
spectrum^ dn/dE « exp(-E/E Q ), where E Q = 150 keV, were used then 
<E> = 94.5 keV. This value is only slightly different from that 
used in generating Table IV. 

It has been demonstrated that energy considerations 
connected with the role of the passive absorber shielding each 
detector are significant if locations are to be determined to 
within a few degrees. A major objective of the BATSE is to make 
those determinations to within less than 1°. It would therefore 
seem necessary to incorporate the effects discussed in the 
present work into the overall scheme of location determinations. 

The difficulties involved in obtaining angular ( 0 , <f> ) 
solutions from a set of statements of the form given by Egn. (7) 
are apparent. The efforts of the present study were restricted 
to obtaining first order corrections to the locations given by 
Eqns. (24) and (25). It has been found possible to obtain the 
second order correction corresponding to the retention of an 
additional term in the expansion of the aluminum exponential 
expansion. The details of that calculation have not yet been 
fully explored, but can and should be done. 
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The need to impose the upper limit on the energy range as 
discussed here to 150 keV resulted from the detector efficiency 
decreasing too rapidly above this energy (as well as somewhat 
below it). Further investigation is clearly necessary so that 
the upper limit can approach the desired 200 keV value. It is 
not entirely clear at this point that those calculations can be 
performed analytically. 
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GAMMA RAY OBSERVATORY 


+X 



BATSE DETECTOR ASPECT 


Figure 1. 


Gamma Ray Observatory Configuration; BATSE Octohedral Geometry 





TABLE I 


Attenuation Coefficients for Aluminum 
and Evaluation of Their Influence 


* 

1 


* 


(keV) 

/ -1\ 
y (cm ) 

e” y 

1-yt 

r 

fi 

50 

0.972 

0.823 

0.806 

k 

75 

0.621 

0.883 

0.876 

i 

100 

0.473 

0.910 

0.905 

k 

125 

0.405 

0.922 

0.919 

it 

150 

0.378 

0.927 

0.924 

f*. 

175 

0.351 

0.932 

0.930 


200 

0.324 

0.937 

0.935 



TABLE II 


Attenuation Coefficients for Sodium Iodide 
and Evaluation of Their Influence 


E(keV) y ( cm" 1 ) e~ yD (1-e pD ) (1-e mD/cos 


50 

36.70 

5.7 x 10~ 21 

>1 

-1 

75 

12.85 

2.6 x 10" 6 

>1 

'1 

100 

6.97 

1.4 x 10" 4 

>1 

'1 

125 

3.67 

9.8 x 10" 3 

0.99 

-1 

150 

2.20 

6.1 x 10“ 2 

0.94 

0.98 

175 

1.47 

0.155 

0.85 

0.93 

200 

1.28 

0.278 

0.72 

0.90 


45° 


) 


XL- 20 


TABLE III 


Computed Angular Deviations of Source as a 
Function of Energy and Fixed N 3 


E(keV) 

T o A 

♦ 

9 

* ( °> 

q(°) 

6 <J> 

60 

50 

4.57 

63.58° 

77.44° 

56.31° 

74.50° 

7.27° 

2.94 

75 

3.96 

60.73° 

76.26° 

56.31° 

74.50° 

4.42° 

1.76 

100 

3.75 

59.61° 

75.80° 

56.31° 

74.50° 

3.30° 

1.30 

125 

3.67 

59.11° 

75.60° 

56.31° 

74.50° 

2.80° 

1.10 

150 

3.64 

58.92“ 

75.53° 

56.31° 

74.50° 

2.61° 

1.03 












TABLE 

IV 





Computed 
for <E> 

Angular Deviations of 
and Various N^, N 2 , N 3 

Source 
( xl 0 n ) 


N 1 

n 2 

n 3 

: o A 

6 $ 

60 

(69 + ) 

3 

2 

1 

3.80 

+3.59° 

+1.42° 

( + ) 

1 

3 

2 

5.68 

-1.31° 

-2.27° 

(-) 

2 

1 

3 

4.42 

-1.57° 

+1.45° 

( + ) 

2 

2 

3 

5.35 

-1.31° 

0 ° 


3 

2 

2 

4.42 

+1.57° 

+1.45° 

( + ) 

2 

3 

2 

5.35 

0 ° 

-1.25° 

(-) 
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